CYBB-Mediated Ferroptosis Associated with Immunosuppression in Mycobacterium leprae-Infected Monocyte-Derived Macrophages.

Mycobacterium leprae-infected macrophages preferentially exhibit the regulatory M2 phenotype in vitro, which helps the immune escape unabated growth of M leprae in host cells. The mechanism that triggers macrophage polarization is still unknown. In this study, we performed single-cell RNA sequencing to determine the initial responses of human monocyte-derived macrophages against M leprae infection of 4 healthy individuals and found an increase in a major alternative-activated macrophage type that overexpressed NEAT1, CCL2, and CD163. Importantly, further functional analysis showed that ferroptosis was positively correlated with M2 polarization of macrophages, and in vitro experiments have shown that inhibition of ferroptosis promotes the survival of M leprae within macrophages. In addition, further joint analysis of our results with mutisequencing data from patients with leprosy and in vitro validation identified that CYBB was the pivotal molecule for ferroptosis that could promote the M2 polarization of M leprae-infected macrophages, resulting in the immune escape and unabated growth of pathogenic bacteria. Overall, our results suggest that M leprae facilitated its survival by inducing CYBB-mediated macrophage ferroptosis leading to its alternative activation and might reveal the potential for a new therapeutic strategy of leprosy.

Genome-wide meta-analysis and fine-mapping prioritize potential causal variants and genes related to leprosy.

To date, genome-wide association studies (GWASs) have discovered 35 susceptible loci of leprosy; however, the cumulative effects of these loci can only partially explain the overall risk of leprosy, and the causal variants and genes within these loci remain unknown. Here, we conducted out new GWASs in two independent cohorts of 5007 cases and 4579 controls and then a meta-analysis in these newly generated and multiple previously published (2277 cases and 3159 controls) datasets were performed. Three novel and 15 previously reported risk loci were identified from these datasets, increasing the known leprosy risk loci of explained genetic heritability from 23.0 to 38.5%. A comprehensive fine-mapping analysis was conducted, and 19 causal variants and 14 causal genes were identified. Specifically, manual checking of epigenomic information from the Epimap database revealed that the causal variants were mainly located within the immune-relevant or immune-specific regulatory elements. Furthermore, by using gene-set, tissue, and cell-type enrichment analyses, we highlighted the key roles of immune-related tissues and cells and implicated the PD-1 signaling pathways in the pathogenetic mechanism of leprosy. Collectively, our study identified candidate causal variants and elucidated the potential regulatory and coding mechanisms for genes associated with leprosy.

Overview and Current Advances in Dapsone Hypersensitivity Syndrome.

PURPOSE OF REVIEW: As a sulfone antibacterial agent, dapsone has been widely used to treat leprosy. Moreover, dapsone is also used in many immune diseases such as herpetic dermatitis because of its anti-inflammatory and immunomodulatory effects. However, dapsone can cause several adverse effects, the most serious being dapsone hypersensitivity syndrome. Dapsone hypersensitivity syndrome is characterized by a triad of eruptions, fever, and organ involvement, which limits the application of dapsone to some extent. RECENT FINDINGS: In this article, we review current research about the interaction model between HLA-B*13:01, dapsone, and specific TCR in dapsone-induced drug hypersensitivity. In addition to the proposed mechanisms, we also discussed clinical features, treatment progress, prevalence, and prevention of dapsone hypersensitivity syndrome. These studies reveal the pathogenesis, clinical features, and prevalence from the perspectives of genetic susceptibility and innate and adaptive immunity in dapsone hypersensitivity syndrome, thereby guiding clinicians on how to diagnose, prevent, and treat dapsone hypersensitivity syndrome.

Single-cell sequencing analysis reveals development and differentiation trajectory of Schwann cells manipulated by M. leprae.

BACKGROUND: M. leprae preferentially infects Schwann cells (SCs) in the peripheral nerves leading to nerve damage and irreversible disability. Knowledge of how M. leprae infects and interacts with host SCs is essential for understanding mechanisms of nerve damage and revealing potential new therapeutic strategies. METHODOLOGY/PRINCIPAL FINDINGS: We performed a time-course single-cell sequencing analysis of SCs infected with M. leprae at different time points, further analyzed the heterogeneity of SCs, subpopulations associated with M. leprae infection, developmental trajectory of SCs and validated by Western blot or flow cytometry. Different subpopulations of SCs exhibiting distinct genetic features and functional enrichments were present. We observed two subpopulations associated with M. leprae infection, a stem cell-like cell subpopulation increased significantly at 24 h but declined by 72 h after M. leprae infection, and an adipocyte-like cell subpopulation, emerged at 72 h post-infection. The results were validated and confirmed that a stem cell-like cell subpopulation was in the early stage of differentiation and could differentiate into an adipocyte-like cell subpopulation. CONCLUSIONS/SIGNIFICANCE: Our results present a systematic time-course analysis of SC heterogeneity after infection by M. leprae at single-cell resolution, provide valuable information to understand the critical biological processes underlying reprogramming and lipid metabolism during M. leprae infection of SCs, and increase understanding of the disease-causing mechanisms at play in leprosy patients as well as revealing potential new therapeutic strategies.

IL-23/IL-23R Promote Macrophage Pyroptosis and T Helper 1/T Helper 17 Cell Differentiation in Mycobacterial Infection.

Pathogen-induced epigenetic modifications can reshape anti-infection immune processes and control the magnitude of host responses. DNA methylation profiling has identified crucial aberrant methylation changes associated with diseases, thus providing biological insights into the roles of epigenetic factors in mycobacterial infection. In this study, we performed a genome-wide methylation analysis of skin biopsies from patients with leprosy and healthy controls. T helper 17 differentiation pathway was found to be significantly associated with leprosy through functional enrichment analysis. As a key gene in this pathway, IL-23R was found to be critical to mycobacterial immunity in leprosy, according to integrated analysis with DNA methylation, RNA sequencing, and GWASs. Functional analysis revealed that IL-23/IL-23R-enhanced bacterial clearance by activating caspase-1/GSDMD-mediated pyroptosis in a manner dependent on NLRP3 through signal transducer and activator of transcription 3 signaling in macrophages. Moreover, IL23/IL-23R promoted T helper 1 and T helper 17 cell differentiation and proinflammatory cytokine secretion, thereby increasing host bactericidal activity. IL-23R knockout attenuated the effects and increased susceptibility to mycobacterial infection mentioned earlier. These findings illustrate the biological functions of IL-23/IL-23R in modulating intracellular bacterial clearance in macrophages and further support their regulatory effects in T helper cell differentiation. Our study highlights that IL-23/IL-23R might serve as potential targets for the prevention and treatment of leprosy and other mycobacterial infections.

Risk factors for dapsone resistance in leprosy patients: a systematic meta-analysis.

OBJECTIVES: Dapsone is one of the important drugs in the treatment of leprosy. The present study aims to evaluate the resistance of Mycobacterium leprae isolates to dapsone, in turn assisting in implementing better control strategies for leprosy elimination. METHODS: A systematic literature search was conducted in PubMed, Embase, Medline, and Web of Science. Two independent reviewers selected the literature according to the guidelines of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA), extracted data, and evaluated the risk of bias. Drug resistance data were pooled using the random-effects model. Subgroup analysis was performed based on across sampling time, region, study population (treatment status, relapses status), and sample size. RESULTS: A total of 30 studies were included. The results of meta-analysis showed that the dapsone resistance rate of leprosy patients after treatment was 8% (95% confidence interval [CI], 6%-10%). Compared to the rates of primary resistance of new cases without treatment therapy (pooled incidence, 4% [95% CI, 2%-5%]), treatment cases (13% [95% CI 9%-16%]) had secondary resistance, and relapse cases (26% [95% CI, 18%-33%]) had drug resistance. In addition, the drug resistance rate of monotherapy was significantly increased than that of relapsed patients treated with diamino-diphenylsulfone monotherapy. Subgroup analysis showed that the patients in the Western Pacific have the highest dapsone resistance, and the resistance to dapsone was slightly lower after 2005. For sample size, the rate in the group under 100 samples was significantly higher than in the other. CONCLUSION: Dapsone resistance is closely related to leprosy relapse and long-term drug use. Dapsone monotherapy is one of important reasons for drug resistance in relapsed cases. Drug resistance varies among different populations and regions of the world.

Utility of Multi-target Nested PCR and ELISPOT Assays for the Detection of Paucibacillary Leprosy: A Possible Conclusion of Clinical Laboratory Misdiagnosis.

The diagnosis of paucibacillary (PB) leprosy often possesses a diagnostic challenge, especially for pure neuritic and lesser skin lesions with the zero bacillary load, requiring a sensitive and accurate diagnostic tool. We have included 300 clinically diagnosed new leprosy cases (comprising 98 PB cases) and analyzed the sensitivity and specificity of PB leprosy cases by nested PCR with folP, gyrA, rpoB, RLEP, and 16SrRNA and Enzyme-linked Immunospot Assay test (ELISPOT) with MMPII, NDO-BSA, and LID-1 antigens by detecting interferon gamma (IFN-γ) release. The overall positivity rates of genes tested in 300 clinical specimens were identified as 55% of 16SrRNA, 59% of RLEP, 59.3% of folP, 57.3% of rpoB, 61% of gyrA while 90% of nested folP, 92.6% of nested rpoB, and 95% of nested gyrA, and 285 (95%) of at least one gene positive cases. For PB specimens, 95% PCR positivity was achieved by three tested genes in nested PCR. The data obtained from ELISPOT for three antigens were analyzed for IFN-γ expression with 600 subjects. Among 98 PB leprosy cases, the sensitivity of MMP II, LID-1, and NDO-BSA was 90%, 87%, and 83%, respectively, and the specificity was 90%, 91%, and 86%, respectively. The total number of cases positive for at least one antigen was 90 (91.8%) in PB, which is significantly higher than that in multibacillary (MB) leprosy (56.7%). The combination of multi-targets nested PCR and ELISPOT assay provides a specific tool to early clinical laboratory diagnosis of PB leprosy cases. The two assays are complementary to each other and beneficial for screening PB patients.

The immune-suppressive landscape in lepromatous leprosy revealed by single-cell RNA sequencing.

Lepromatous leprosy (L-LEP), caused by the massive proliferation of Mycobacterium leprae primarily in macrophages, is an ideal disease model for investigating the molecular mechanism of intracellular bacteria evading or modulating host immune response. Here, we performed single-cell RNA sequencing of both skin biopsies and peripheral blood mononuclear cells (PBMCs) of L-LEP patients and healthy controls. In L-LEP lesions, we revealed remarkable upregulation of APOE expression that showed a negative correlation with the major histocompatibility complex II gene HLA-DQB2 and MIF, which encodes a pro-inflammatory and anti-microbial cytokine, in the subset of macrophages exhibiting a high expression level of LIPA. The exhaustion of CD8+ T cells featured by the high expression of TIGIT and LAG3 in L-LEP lesions was demonstrated. Moreover, remarkable enhancement of inhibitory immune receptors mediated crosstalk between skin immune cells was observed in L-LEP lesions. For PBMCs, a high expression level of APOE in the HLA-DRhighFBP1high monocyte subset and the expansion of regulatory T cells were found to be associated with L-LEP. These findings revealed the primary suppressive landscape in the L-LEP patients, providing potential targets for the intervention of intracellular bacteria caused persistent infections.

Single-Nucleotide Polymorphisms Related to Leprosy Risk and Clinical Phenotypes Among Chinese Population.

BACKGROUND: Genome-wide association studies (GWASs) have identified some immune-related single-nucleotide polymorphisms (SNPs) to be associated with leprosy. METHODS: This study investigated the association of 17 SNPs based on previously published GWAS studies with susceptibility to leprosy, different polar forms and immune states of leprosy in a case-control study from southwestern China, including 1344 leprosy patients and 2732 household contacts (HHCs) (1908 relatives and 824 genetically unrelated contact individuals). The differences of allele distributions were analyzed using chi-squared analysis and logistic regression. RESULTS: After adjusting covariate factors, rs780668 and rs3764147 polymorphisms influenced susceptibilities to genetically related or unrelated leprosy contact individuals. rs142179458 was associated with onset early cases, rs73058713 A allele and rs3764147 A allele increased the risk of reversal reaction, while rs3764147 G allele had higher risk to present lepromatous leprosy and erythema nodosum leprosum. CONCLUSION: Our results demonstrated that genetic variants in the LACC1, HIF1A, SLC29A3 and CDH18 genes were positively correlated with the occurrence of leprosy and leprosy clinical phenotypes, providing new insights into the immunogenetics of the disease.

Predictive nomogram for leprosy using genetic and epidemiological risk factors in Southwestern China: Case-control and prospective analyses.

BACKGROUND: There is a high incidence of leprosy among house-contacts compared with the general population. We aimed to establish a predictive model using these genetic factors along with epidemiological factors to predict leprosy risk of leprosy household contacts (HHCs). METHODS: Weighted genetic risk score (wGRS) encompassing genome wide association studies (GWAS) variants and five non-genetic factors were examined in a case-control design associated with leprosy risk including 589 cases and 647 controls from leprosy HHCs. We constructed a risk prediction nomogram and evaluated its performance by concordance index (C-index) and calibration curve. The results were validated using bootstrap resampling with 1000 resamples and a prospective design including 1100 HHCs of leprosy patients. FINDING: The C-index for the risk model was 0·792 (95% confidence interval [CI] 0·768-0·817), and was confirmed to be 0·780 through bootstrapping validation. The calibration curve for the probability of leprosy showed good agreement between the prediction of the nomogram and actual observation. HHCs were then divided into the low-risk group (nomogram score ≤ 81) and the high-risk group (nomogram score > 81). In prospective analysis, 12 of 1100 participants had leprosy during 63 months' follow-up. We generated the nomogram for leprosy in the validation cohort (C-index 0·773 [95%CI 0·658-0·888], sensitivity75·0%, specificity 66·8%). Interpretation The nomogram achieved an effective prediction of leprosy in HHCs. Using the model, the risk of an individual contact developing leprosy can be determined, which can lead to a rational preventive choice for tracing higher-risk leprosy contacts. FUNDING: The ministry of health of China, ministry of science and technology of China, Chinese academy of medical sciences, Jiangsu provincial department of science and technology, Nanjing municipal science and technology bureau.

Massively Parallel Sequencingof the Filaggrin Gene Reveals an Association Between FLG Loss-of-function Mutations and Leprosy.

Filaggrin, encoded by the FLG gene, plays a crucial role in the barrier function of epidermis, but the association between FLG loss-of-function mutations and infectious skin diseases has not been systematically studied. FLG coding sequences from 945 patients with leprosy and 916 healthy controls were captured and enriched using an array-based high-throughput system, and subjected to next-generation sequencing. The loss-of-function mutations found were further validated by Sanger sequencing. A total of 21 loss-of-function mutations were found in 945 patients with leprosy, with a carrier rate of 17.53%, while the prevalence of these mutations in 916 healthy controls was 14.77%, which was significantly lower than in patients. Two individual FLG loss-of-function mutations (K4022X and Q1790X) were found to be significantly associated with leprosy. These results suggest a possible role for filaggrin in defending against leprosy pathogens.

The role of an active surveillance strategy of targeting household and neighborhood contacts related to leprosy cases released from treatment in a low-endemic area of China.

OBJECTIVE: Early diagnosis remains the primary goal for leprosy management programs. This study aims to determine whether active surveillance of patients with leprosy and their contact individuals increased identification of latent leprosy cases in the low-endemic areas. METHODS: This cross-sectional survey was carried out between October 2014 and August 2016 in 21 counties throughout Shandong Province. The survey was conducted among patients with leprosy released from treatment (RFT) and their contacts from both household and neighbors. RESULTS: A total of 2,210 RFT patients and 9,742 contacts comprising 7877 household contacts (HHCs), including 5,844 genetic related family members (GRFMs) and 2033 non-genetic related family members and 1,865 contacts living in neighboring houses (neighbor contacts, NCs), were recruited. Among identified individuals, one relapsed and 13 were newly diagnosed, giving a detection rate of 0.12%, corresponding to 120 times the passive case detection rate. Detection rates were similar for HHCs and NCs (0.114% vs. 0.214%, P = 0.287). Analysis of the family history of leprosy patients revealed clustering of newly diagnosed cases and association with residential coordinates of previously-diagnosed multibacillary leprosy cases. CONCLUSION: Active case-finding programs are feasible and contributes to early case detection by tracking HHCs and NCs in low-endemic areas.

Evaluation of Prospective HLA-B*13:01 Screening to Prevent Dapsone Hypersensitivity Syndrome in Patients With Leprosy.

Importance: Dapsone hypersensitivity syndrome (DHS) is the most serious adverse reaction associated with dapsone administration and one of the major causes of death in patients with leprosy, whose standard treatment includes multidrug therapy (MDT) with dapsone, rifampicin, and clofazimine. Although the HLA-B*13:01 polymorphism has been identified as the genetic determinant of DHS in the Chinese population, no studies to date have been done to evaluate whether prospective HLA-B*13:01 screening could prevent DHS by identifying patients who should not receive dapsone. Objective: To evaluate the clinical use of prospective HLA-B*13:01 screening for reduction of the incidence of DHS by excluding dapsone from the treatment for patients with HLA-B*13:01-positive leprosy. Design, Setting, and Participants: A prospective cohort study was conducted from February 15, 2015, to April 30, 2018, in 21 provinces throughout China. A total of 1539 patients with newly diagnosed leprosy were enrolled who had not received dapsone previously. After excluding patients who had a history of allergy to sulfones or glucose-6-phosphate dehydrogenase deficiency, 1512 individuals underwent HLA-B*13:01 genotyping. All of the patients were followed up weekly for the first 8 weeks after treatment to monitor for adverse events. Exposures: Patients who were HLA-B*13:01 carriers were instructed to eliminate dapsone from their treatment regimens, and noncarrier patients received standard MDT. Main Outcomes and Measures: The primary outcome was the incidence of DHS. The historical incidence rate of DHS (1.0%) was used as a control. Results: Among 1512 patients (1026 [67.9%] men, 486 [32.1%] women; mean [SD] age, 43.1 [16.2] years), 261 (17.3%) were identified as carriers of the HLA-B*13:01 allele. A total of 714 adverse events in 384 patients were observed during the follow-up period. Dapsone hypersensitivity syndrome did not develop in any of the 1251 patients who were HLA-B*13:01-negative who received dapsone, while approximately 13 patients would be expected to experience DHS, based on the historical incidence rate of 1.0% per year (P = 2.05 × 10-5). No significant correlation was found between other adverse events, including dermatologic or other events, and HLA-B*13:01 status. Conclusions and Relevance: Prospective HLA-B*13:01 screening and subsequent elimination of dapsone from MDT for patients with HLA-B*13:01-positive leprosy may significantly reduce the incidence of DHS in the Chinese population.

Development and evaluation of a droplet digital PCR assay for the diagnosis of paucibacillary leprosy in skin biopsy specimens.

BACKGROUND: The reduced amounts of Mycobacterium leprae (M. leprae) among paucibacillary (PB) patients reflect the need to further optimize methods for leprosy diagnosis. An increasing number of reports have shown that droplet digital polymerase chain reaction (ddPCR) is a promising tool for diagnosis of infectious disease among samples with low copy number. To date, no publications have investigated the utility of ddPCR in the detection of M. leprae. The aim of this study was to develop and evaluate a ddPCR assay for the diagnosis of PB leprosy. METHODOLOGY: The two most sensitive DNA targets for detection of M. leprae were selected from electronic databases for assessment of sensitivity and specificity by quantitative polymerase chain reaction (qPCR) and ddPCR. Control patients (n = 59) suffering from other dermatological diseases were used to define the cut-off of the duplex ddPCR assay. For comparative evaluation, qPCR and ddPCR assays were performed in 44 PB patients and 68 multibacillary (MB) patients. PRINCIPAL FINDINGS: M. leprae-specific repetitive element (RLEP) and groEL (encoding the 65 kDa molecular chaperone GroEL) were used to develop the ddPCR assay by systematically analyzing specificity and sensitivity. Based on the defined cut-off value, the ddPCR assay showed greater sensitivity in detecting M. leprae DNA in PB patients compared with qPCR (79.5% vs 36.4%), while both assays have a 100% sensitivity in MB patients. CONCLUSIONS/SIGNIFICANCE: We developed and evaluated a duplex ddPCR assay for leprosy diagnosis in skin biopsy samples from leprosy patients. While still costly, ddPCR might be a promising diagnostic tool for detection of PB leprosy.

Prediction of leprosy in the Chinese population based on a weighted genetic risk score.

Genome wide association studies (GWASs) have revealed multiple genetic variants associated with leprosy in the Chinese population. The aim of our study was to utilize the genetic variants to construct a risk prediction model through a weighted genetic risk score (GRS) in a Chinese set and to further assess the performance of the model in identifying higher-risk contact individuals in an independent set. The highest prediction accuracy, with an area under the curve (AUC) of 0.743 (95% confidence interval (CI): 0.729-0.757), was achieved with a GRS encompassing 25 GWAS variants in a discovery set that included 2,144 people affected by leprosy and 2,671 controls. Individuals in the high-risk group, based on genetic factors (GRS > 28.06), have a 24.65 higher odds ratio (OR) for developing leprosy relative to those in the low-risk group (GRS≤18.17). The model was then applied to a validation set consisting of 1,385 people affected by leprosy and 7,541 individuals in contact with leprosy, which yielded a discriminatory ability with an AUC of 0.707 (95% CI: 0.691-0.723). When a GRS cut-off value of 22.38 was selected with the optimal sensitivity and specificity, it was found that 39.31% of high risk contact individuals should be screened in order to detect leprosy in 64.9% of those people affected by leprosy. In summary, we developed and validated a risk model for the prediction of leprosy that showed good discrimination capabilities, which may help physicians in the identification of patients coming into contact with leprosy and are at a higher-risk of developing this condition.

Amino Acid Variants of HLA-DRB1 Confer Susceptibility to Dapsone Hypersensitivity Syndrome in Addition to HLA-B*13:01.

Dapsone hypersensitivity syndrome is a rare yet severe adverse drug reaction caused by dapsone, a principal drug in multidrug therapy for leprosy. HLA-B*13:01 has been identified as a strong risk factor of dapsone hypersensitivity syndrome; however, its low positive predictive value indicated that additional genetic variants may be involved in the disease development. To discover contributing genetic variants within HLA loci in addition to HLA-B*13:01, we performed a high-coverage next-generation sequencing (NGS)-based HLA typing analysis in 103 dapsone-hypersensitive and 857 dapsone-tolerant HLA-B*13:01-positive leprosy patients in a Chinese population. Five amino acid variants in high linkage disequilibrium of HLA-DRB1 were significantly associated with dapsone hypersensitivity syndrome (positions 133, 142, -17, 11, and 13). DRB1*16:02 and DRB1*15:01 tagged by these risk-conferring amino acid residues were associated at a nominal significance level. This study identifies five amino acid variants within HLA-DRB1 that are in high linkage disequilibrium and significantly associated with dapsone hypersensitivity syndrome in a Chinese population.

A pathway-based association analysis identified FMNL1-MAP3K14 as susceptibility genes for leprosy.

The nuclear transcription factor-κB (NF-κB) plays a pivotal role in controlling both innate and adaptive immunity and regulates the expressions of many immunological mediators. Abundant evidences have showed the importance of NF-κB pathway in the host immune responses against Mycobacterium leprae in the development of leprosy. However, no particular association study between leprosy and NF-κB pathway-related gene polymorphisms was reported. Here, we performed a large-scale and two-stage candidate association study to investigate the association between 94 NF-κB pathway-related genes and leprosy. Our results showed that rs58744688 was significantly associated with leprosy (P = 7.57 × 10-7 , OR = 1.12) by combining the previous genomewide association data sets and four independent validation sample series, consisting of a total of 4631 leprosy cases and 6413 healthy controls. This founding implicated that MAP3K14 and FMNL1 were susceptibility genes for leprosy, which suggested the involvement of macrophage targeting and NF-κB pathway in the development of leprosy.

Genome-Wide Analysis of Protein-Coding Variants in Leprosy.

Although genome-wide association studies have greatly advanced our understanding of the contribution of common noncoding variants to leprosy susceptibility, protein-coding variants have not been systematically investigated. We carried out a three-stage genome-wide association study of protein-coding variants in Han Chinese, of whom were 7,048 leprosy patients and 14,398 were healthy control subjects. Seven coding variants of exome-wide significance were discovered, including two rare variants: rs145562243 in NCKIPSD (P = 1.71 × 10-9, odds ratio [OR] = 4.35) and rs149308743 in CARD9 (P = 2.09 × 10-8, OR = 4.75); three low-frequency variants: rs76418789 in IL23R (P = 1.03 × 10-10, OR = 1.36), rs146466242 in FLG (P = 3.39 × 10-12, OR = 1.45), and rs55882956 in TYK2 (P = 1.04 × 10-6, OR = 1.30); and two common variants: rs780668 in SLC29A3 (P = 2.17 × 10-9, OR = 1.14) and rs181206 in IL27 (P = 1.08 × 10-7, OR = 0.83). Discovered protein-coding variants, particularly low-frequency and rare ones, showed involvement of skin barrier and endocytosis/phagocytosis/autophagy, in addition to known innate and adaptive immunity, in the pathogenesis of leprosy, highlighting the merits of protein-coding variant studies for complex diseases.